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The additivity principle allows a calculation of current fluctuations and associated density profiles 
in large diffusive systems. In order to test its validity in the weakly asymmetric exclusion process 
with open boundaries, we use a numerical approach based on the density matrix renormalisation. 
With this technique, we determine the cumulant generating function of the current and the density 
profile corresponding to atypical currents in finite systems. We find that these converge to those 
predicted by the additivity principle. No evidence for dynamical phase transitions is found. 
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Systems driven out of equilibrium by putting them 
in contact with reservoirs at different chemical poten- 
tial or temperature develop currents. It is a main prob- 
lem of non equilibrium statistical mechanics to deter- 
mine how the average and fluctuations of these currents 
can be derived from microscopic dynamics. In recent 
years, considerable progress has been made in this di- 
rection. Firstly, it was found that current fluctuations 
have symmetries as expressed in the Gallavotti-Cohen 
[TJ [2] theorem. These symmetries are macroscopic man- 
ifestations of microscopic time reversibility. Secondly, 
Bodineau and Derrida formulated an additivity princi- 
ple (AP) that allows one to calculate the whole distribu- 
tion of current fluctuations once the first two cumulants 
are known j3]. The AP should hold for one-dimensional, 
diffusive systems and was validated in the symmetric ex- 
clusion processes [I] and the Kipnis-Marchioro-Presutti 
(KMP) model of heat conduction [5]. Most recently, it 
was also found to hold in three-dimensional determinis- 
tic models of heat conduction [6 . Independently, Bertini 
et al. developped a large deviation theory for density 
and current fluctuations in stochastic lattice gases [7, !8J. 
The predictions of this Hydrodynamic Fluctuation The- 
ory coincide with those of the AP when the fluctuations 
are time-independent. Interestingly, it was found that for 
sufficiently large fluctuations a dynamical phase transi- 
tion can occur to a phase where density and current fluc- 
tuations become time-dependent [5] . This transition was 
observed in a weakly asymmetric exclusion process on a 
ring [9] and more recently in the KMP model [10], also 
on a ring. This type of dynamical transitions can oc- 
cur in situations that are not allowed in equilibrium and 
have been conjectured to be of relevance to such issues 
as breaking of chiral or CP-symmetry [TT] . 

In this Letter, we study the current fluctuations in 
the weakly asymmetric exclusion process (WASEP) with 
open boundaries. We calculate the cumulant generat- 
ing function of the current and the density profile giving 
rise to an atypical current using the AP. This extends 
earlier work [IS]. We compare these results with those 



coming from calculations in finite systems using the den- 
sity matrix renormalisation group (DMRG). We recently 
showed how that approach, first introduced to study low 
temperature properties of quantum systems [13], can be 
applied to determine the cumulant generating function 
of the current or the activity of stochastic systems [14] . 
In the present work, we are the first to calculate den- 
sity profiles corresponding to current fluctuations with 
the DMRG. We find that for sufficiently large systems 
our results converge to those predicted by the AP in the 
whole range of parameters investigated, further validat- 
ing this principle. We find no evidence for a dynamical 
transition in this case. 

Model — In the asymmetric exclusion process (ASEP) 
[4] , each site i of a lattice of size N can be empty or occu- 
pied by one particle. The dynamics is that of a Markov 
chain where particles jump to the right or left with dif- 
ferent rates (p resp. q). This describes the effect of an 
external field E = ln(p/q). In this Letter we will discuss 
the case where p = 1 + v/ (2N), q = 1 — vj (2N) referred 
to as the WASEP. This is a diffusive model where the 
AP should be applicable. We will consider the case of 
open boundaries where at its left (right) side the system 
is in contact with a reservoir at density p a (pb)- We will 
assume p a > pb and v > 0. The model will evolve to a 
non equilibrium steady state (NESS). 

Current fluctuations from the additivity principle — We 
are interested in the total number of particles Qt passing 
in a large time T through the system. For a large system, 
and using a continuum description in terms of x = i/N € 
[0, 1] , the average current equals 

(Q T ) = D(p) dp vo-(p) f 
T N dx N N 1 ' 

The first term is Fick's law and the second is the current 
due to the field E in linear response. For the WASEP, the 
diffusivity D(p) = 1/2 and the mobility a(p) = p(l — p) 
[15j where p(x) is the particle density. In this Letter, we 
are interested in the fluctuations of the current around 
the average value j*. For N very large, the probability 
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Pn to observe an integrated current Qt = jT/N has the 
form 



Qt 
T 
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G(j,Pa,Pb) 



(2) 



The large deviation function G is zero at the average 
current j* of the NESS, and is strictly negative for other 
j-valucs. According to the AP [3], G can be found from 
a variational principle 



G(j,Pa,Pb) 
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>o 1a(p{x)) 
which leads to a Euler-Lagrange equation for p(x) 
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= (j - va{p)f + 2Ka(p) 



(4) 



Solution of this equation for given j and subject to the 
boundary conditions p(0) = p a , p(l) = Pb determines 
the integration constant K and the density profile. For 
a monotonically decreasing profile, one obtains 



dp 



Pa [{ J ~ V a{pW + 2Ko{p)] 1 l 2 



= -2x 



(5) 



More complicated expressions can be determined for the 
case that the profile has an extremum. Inserting this 
solution in ([3| then gives for the large deviation function 



G(j,p a ,Pb) = 



(j - va{p)) 



(j - va{p)f + Ka(p) 
[{j-ya{p)Y + 2Ka{p)]^ 
dp 



Mp) 



(6) 



Instead of ([2|, one can also describe the current fluc- 
tuations using the cumulant generating function 

Ks,Pa,Pb)= lim ^ln(e s< 3-) (7) 
which is related to G through a Legendre transform 

i \ 1 r • . nc w M ( s ?Pa,Pb) , Q , 

P{s, Pa, Pb) = ^7 max [sj + G(j, p a , p b )\ = (8) 

Inserting ^ gives the cumulant generating function 
(CGF) in parametric form 



M(s,p a ,p b ) 



1_ /""• (j - va{p))v ~K d 
2./ Pb [{ ] -ua{p)Y+2Ka{p)] 1/2 '' 



v 



and 



'Pb 



+ ^{Pb-Pa) 



j - va(p) 
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This is again the result for a monotonically decreasing 
profile. The more complicated expressions for a profile 
with an extremum will be given elsewhere. For given 
values of is, p a , pt, and j we have determined the density 
profile, the large deviation function and the CGF by nu- 
merical evalution of the integrals in ([5]), ^ and (10 1. 

DMRG - approach — We now want to determine the 
density profile and the CGF for finite systems to see 
whether for large N they converge to those predicted by 
the AP. With standard simulation techniques it is diffi- 
cult to generate atypical currents since they occur with 
exponentially small probability. A method to overcome 
this problem has been proposed in [TU [17]. Yet, this 
technique becomes less accurate for large fluctuations [5] 
due to statistical errors. Recently, we proposed a new 
approach to current fluctuations based on the DMRG, 
the ideas behind which we now briefly explain |18| . 

The probability P(C,t) to observe the exclu- 
sion process in a given microscopic configuration C 
evolves according to the master equation d t P(C,t) = 
^2, c , H(C, C')P(C, t) where H is the generator of the pro- 
cess. It is by now well established g] that the CGF Q 
can be obtained from a modified generator H s , which 
is constructed from H as follows. Let a(C,C) be +1 
( — 1) when in the transition from C to C a particle 
enters (leaves) the system on its left side. Otherwise 
a(C,C) — 0. For the off-diagonal elements of H s one has 
H S (C,C) = iJ(C,C')e SQ(C ' C,) while the diagonal elements 
of H and H s are equal. The CGF for a system of N sites, 
p(s, p a , pb, N), then equals the largest eigenvalues of H s . 
Moreover, let |i?o) an d (^ol be the associated right and 
left eigenvector. Consider a dynamical variable b(C(t)) 
(like the density at a given site) which depends on the 
microscopic configuration in which the system is at time 
t. It can be shown that the current weighted time-average 
of b, defined as, 



(hM) _l (fib(C(T))e^ dT) 



(11) 



for T large equals (L \b\R ) where b is the operator as- 
sociated to the variable b [TS]. On the other hand, the 
average of b at a large time T defined as 



(b(T)) s = (&(C(T))e^)/(e 



sQt\ 



(12) 



equals (0\b\R o ) where (0| is the projection state (C| . 
So both the CGF and the density profiles can be de- 
termined from the largest eigenvalue of H(s), its eigen- 
vectors and the projection state. From a mathematical 
point of view, solving this problem is similar to that of 
determining the ground state and its eigenvector for a 
quantum spin chain, where the main difference is that in 
the stochastic problem the generator H s is not Hcrmi- 
tian. One of the most precise approaches to determine 
ground state properties of quantum chains is the DMRG 
[T31 |2"U] . We recently showed that this method also works 
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FIG. 1: (Color online) Cumulant generating for the WASEP 
for v — 10, p a = 4/7, pb = 5/18 from the additivity principle 
and from the DMRG. 



FIG. 2: (Color online) Current large deviation function for 
the WASEP for v = 10, p a = 4/7, p b = 5/18 from the additiv- 
ity principle (full line). The dotted line indicates a quadratic 
fit near the maximum. 



well for generalized generators associated with current 
fluctuations [14]. Here we give for the first time results 
on time-averaged density profiles. Since the projection 
state plays no role in quantum mechanical problems, we 
had to adapt the DMRG approach in order to also cal- 
culate late-time averages. Details of this will be given 
elsewhere. We can typically obtain reliable results up to 
N w 120. As a check of the DMRG approach we have 
calculated density profiles for the totally asymmetric ex- 
clusion process (p — 1, q — 0) at s — where exact results 
for finite N exist [5T] and have found perfect agreement. 
We are not aware of any exact results for the density of 
the ASEP in finite systems and for s/0. 

Results — We have performed most of our calculations 
for v = 10, p a = 4/7 and p = 5/18. In Figjljwe present 
our results for the CGF. The full line is the prediction 
form the AP. In the regime between the vertical dot- 
ted lines, the optimal profile has a minimum, otherwise 
it is monotonically decreasing. The various symbols are 
DMRG results for N[i(s, p a ,Pb, N) at different N- values. 
As can be seen, for increasing system size, the AP and 
DMRG results coincide within numerical accuracy in an 
increasing range of s- values. Within the whole s-region 
investigated our numerical data satisfy the Gallavotti- 
Cohen symmetry p(s, p a , p b , N) = p(A - s,p a ,p bl N) 
where 

A = - k t^-(iV-l)lni^ (13) 



PaO- ~ Pb) 



l + v/2N 



A dynamical phase transition should show up as a point 
where the CGF becomes non-analytical. On the scale 
of Fig. [I] this seems to occur where the profile changes 
from monotonic to one with a minimum. A detailed in- 
vestigation of the first and second derivative of the CGF 
near these points however shows no evidence for non- 
analyticity. 

In Fig. [2] we show the large deviation function G(j) 
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FIG. 3: (Color online) Time-averaged density profile in a 
finite system N = 50 with current j — ±1.5 as calculated 
with the DMRG. 



as calculated from the AP. This quantity cannot be de- 
termined from the DMRG. For small deviations from the 
average current j*(= 2.5845, as follows from ^ and the 
chosen boundary values), we expect P(j/N) to be Gaus- 
sian, so that the LDF is quadratic (dotted line) . Clearly, 
sufficiently large fluctuations are non-Gaussian. 

We now turn to the density profiles. Firstly, we observe 
that the time-averaged density profiles are invariant for 
the transformation j -h- —j. This is a consequence of 
time-reversibility, and is a special case of a more gen- 
eral result for isometric current fluctuations that holds 
in higher dimensions [22] . The invariance of the profiles 
follows directly from the equations of the AP, but is also 
valid for finite systems. Fig. [3]shows a result for a system 
with N = 50 and j = ±1.5. Density profiles calculated 
after a large time T do not obey this symmetry. 

For various values of s (or j) we have calculated den- 
sity profiles in finite systems. We show as an example in 
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FIG. 4: (Color online) Time-averaged density profile at s = 
10 as calculated from the AP (full line) and from the DMRG. 
The inset shows the density at the first site as a function of 
1/N together with the prediction of the AP (full square). 
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FIG. 5: (Color online) Time-averaged density profile at s = 
— 10 as calculated from the AP (full line) and from the DMRG. 
The inset shows a the time-averaged and the late time profile 
for N = 120. 



Fig. [4] the time-averaged profile corresponding to a large 
positive current fluctuation, s = 10 (or j = 5.1214..). 
The full line is the result from the AP, the symbols indi- 
cate DMRG data. As can be seen, the finite size results 
converge again to those predicted by the AP. This con- 
vergence is slowest near the boundaries. We therefore 
show in the inset an extrapolation of the average den- 
sity at the leftmost side (i = 1) for various N, which 
is consistent with the asymptotic prediction. The figure 
also shows that, apart from boundary effects, the den- 
sity profile becomes flat and concentrated near p = 1/2 
in order to carry this current which is almost twice as 
large as the average one. In Fig. [5] we show similarly 
the density profile associated with a very small current, 
s = —10, j = 0.00041... Also for this profile with a mini- 
mum, finite N results converge to the predictions of the 
AP. In this case the density becomes almost zero, except 
near the boundaries. In the inset we compare the time- 
averaged and late T-profile for a system with N = 120. 
Late time profiles cannot be obtained from the AP. They 
are different from the time-averaged ones in the same way 
as spatial boundaries give rise to differences between bulk 
and surface densities [TS]. 

Conclusions — We have shown that current and density 
fluctuations in the WASEP are for N sufficiently large 
precisely given by the AP. In contrast to the case on a ring 
[9] no evidence for a dynamic phase transition was found. 
This can be seen as a sort of non-equivalence between 
'ensembles' since on a ring particle number is fixed ('mi- 
cro canonical') and for open boundaries it is not ('grand 
canonical'). It is well established that in equilibrium sys- 
tems with long range interactions non-equivalence of en- 
sembles can appear |24j . Non-equilibrium systems like 
the WASEP have long range correlations in time [4] and 
may therefore show similar phenomena. 

From a comparison between the asymptotic AP results 



and those from the DMRG it is possible to quantify finite 
size corrections. This can lead to a finite size scaling 
theory along the lines existing for the totally asymmetric 
exclusion process [T3J [53] . 

We have shown that the DMRG can give reliable re- 
sults on density profiles in systems carrying a large fluc- 
tuation. It can therefore be used with confidence in non 
diffusive models or reaction-diffusion systems where so 
far few analytical approaches to large deviations exist. 
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